Developing the Deutsch-Hayden approach to quantum mechanics 



O 

o 

(N 

> 
O 

m 
> 
in 

oo ' 

o ■ 

On : 
O 

o : 
i!^ • 

Oh! 

> ■ 

a ' 



X 



C. Hewitt-Horsmant, V. Vedral* 
^ Blackett Laboratory, Imperial College, London SW7 2BZ, United Kingdom 
^School of Physics and Astronomy, University of Leeds, Leeds LS2 9JT, United Kingdom 

(Dated: 13 September 2006) 

The formalism of Deutsch and Hayden is a useful tool for describing quantum mechanics explicitly 
as local and unitary, and therefore quantum information theory as concerning a "flow" of information 
between systems. In this paper we show that these physical descriptions of flow are unique, and 
develop the approach further to include the measurement interaction and mixed states. We then 
give an analysis of entanglement swapping in this approach, showing that it does not in fact contain 
non-local effects or some form of superluminal signalling. 

PACS numbers: 



I. INTRODUCTION 

The Deutsch-Hayden approach [l| was introduced in 
order to describe the flow of information through quan- 
tum mechanical systems. The standard formalism is not 
suitable in this situation as it has non-local states; that is, 
a single state can describe systems which are physically 
far apart. Any information stored in that state can then 
appear to 'jump' from one system to the other without 
an explicit physical communication happening. Deutsch 
and Hayden demonstrated that this does not happen in 
their approach, and that only local interactions change 
the description of a system. Furthermore, unlike in the 
standard formalism, there is no notion of "collapse" in 
the Deutsch-Hayden approach. All evolution is unitary, 
even under the action of measurement. This approach is 
therefore very useful if we wish a formalism of quantum 
mechanics that is explicitly local and unitary. One par- 
ticular use of such a formalism is in information theory. 
As Deutsch and Hayden showed in their analysis of tele- 
portation, if we have no notion of collapse then 'bits' and 
'classical communication' are treated by the formalism in 
exactly the same way as qubits and communication via 
quantum channels. This is very useful as we no longer 
need to swap between different types of entities in the 
middle of a protocol. 

Deutsch and Hayden demonstrated information flow in 
some situations, for example in teleportation. However, 
there were some situations that they did not address, and 
some issues that need to be resolved. Firstly, the demon- 
stration of information flow depends on the exact form of 
the description of the systems. However, if we are to use 
this then we need to know how unique those forms are: 
if another form may equally well be used, where does 
that leave the analysis? Then there are issues of mea- 
surement. The measurement interaction itself is unitary, 
but eventually in order to extract prediction from the 
theory we are going to need something that will stand 
in for collapse - a notion of what the description of one 
system is relative to a description of the other. Finally, 
Deutsch and Hayden dealt only with pure state systems, 
so in order to use the formalism universally we are going 
to need a way of describing mixed state systems within 



the approach. 

These, then, are the issues that this paper will address. 
We will first give an introduction to Deutsch and Hay- 
den's work. We will then look at the uniqueness of the 
descriptions given. Then we will deal with the measure- 
ment interaction and relative states. After this, we will 
look in detail at how mixed state systems may be incorpo- 
rated into the formalism. We will then use the developed 
formalism to analyse information flow in a protocol that 
could not be described fully previously: that of entangle- 
ment swapping. 



II. THE DEUTSCH-HAYDEN APPROACH 

The Deutsch-Hayden approach is based on the formal- 
ism flrst introduced by Gottesman Q in the context of 
stabiliser theory, and is a Heisenberg-type representation 
of quantum mechanics. Instead of the states of systems 
evolving under a Hamiltonian and operators being fixed, 
there is a fixed universal state (by convention usually 
0)), and an operator A{t) evolves over time with equa- 
tions of motion 



<:\A{t) 
At 



[A{t),H] 



where H is the Hamiltonian. We can also write the time- 
dependent operator as A{t) — U"^ A{fS)\J , where U is the 
unitary evolution operator. 

The Deutsch-Hayden representation in based on the 
Hilbert space where these time-dependent operators are 
vectors rather than second rank tensors - the space 
of Hilbert-Schmidt operators [3]. This Hilbert-Schmidt 
space has an inner product {^xi^J^B) where A and B 
are operators) and a norm (•\/Tr(A^)), and for an A^- 
dimensional system an A^^ dimensional space is required. 

As an example, let us consider a single qubit. Such a 
system is 2 dimensional and will therefore need a set of 
4 basis operators. One such set that would be useful is 
that of the 2 dimensional Pauli operators, {l,cr}. Now, 
the 1 component of any operator can never change (it 
will evolve as C/^1J7 — 1), and in order for the opera- 
tor to be normalised, the 1 component must always be 



2 



Thus a qubit can be characterised by giving the 
components in the (Jr^,ay^z directions only. This is, of 
course, the Bloch sphere [4, pl74]. 

In the general case, we choose a set of basis operators 
written here as F^. Then any operator can be written in 
terms of the basis (here at time t = 0): 



where S IR is the inner product Tr{A{0)T i{0)) (we as- 
sume that all the Hilbert-Schmidt operators with which 
we work are self-adjoint; that is, they may be the op- 
erators corresponding to observables). If we now evolve 
A{0) to time t we have 

A{t) = U^A{Q)U= J2 a^r,(i) 

i=0 

Thus we see that in order to find the time-evolved state of 
an operator, all we need do is find the time-evolved basis 
operators and then reconstruct the operator using the 
original coefficients. Thus a complete characterisation of 
the time-evolved system can be gained by following the 
time evolution of the basis operators. 

However, this is not usually very practical: as N in- 
creases, the task of following the evolution of iV^ basis op- 
erators quickly becomes unwieldy. It was Gottesman Q 
who realised that the number of operators to be tracked 
can be reduced further. Operators can always be writ- 
ten as combinations of other operators, either additively 
or multipicatively. It is obvious that for additive groups 
of operators the evolution operator preserves the group 
structure (le U^A + B)U = WAU + U^BU). However 
it is also true of the multiplicative groups, where U forms 
a group homomorphism between {X} and {C/^XC/}: 

A{Q)B(Q) U^ABU = U'^AUU^BU 
= A{t)B{t) 

We can therefore see that, rather than following the 
complete set of basis operators, all we need follow is a 
generating set of the group. The time-evolved generating 
set will preserve the group structure of the complete set, 
enabling that to be constructed out of it. The generating 
set that Gottesman (whom Deutsch and Hayden follow) 
chooses is the Pauli group. This gives us the representa- 
tion of an n-qubit array, each qubit of which is defined 
by a set of operators that will generate all the operators 
pertaining to that qubit. We note that this representa- 
tion does in fact have a further redundancy within itself, 
as the 2-dimensional Pauli operators themselves form a 
multiplicative group {eg ay = lUxaz)- 

When using this formalism, not only do we need to 
time-evolve the descriptors, qa, we also need to use time- 
evolved forms of transformation operators. These forms 
can be written in terms of their actions on the descriptors 



over a time-step. For example, the Hadamard gate is 
written in its general form as 

q(ii) = I 

qx{t) 

and the CNOT gate (with qubit 1 as control and 2 as 
target) 

q2x{t)qix{t) \ I q2x{t) 

qi(ii)= ( q2x{t)qiy{t) q2 = q2y{t)qiz{t) 
qiz{t) ) \q2z{t)qiz{t) 

Two useful operators are 

z± = \{t±qz) 

Their average values are 

(z+)=Tr(p|0}(0|) (z_)=Tr(Hl)(l|) 

That is, the average values of these operators are the 
diagonal elements of the density operator in the compu- 
tational basis. 

The use of these operators scales very easily as we add 
more systems. We saw above how the qai for each system 
always describe separate subspaces for the individual sys- 
tem. We can therefore simply combine the z± for each 
system to get overall probabilities. For example, if we 
have two systems then 

p(01) = (Zi+®Z2_) 

which gives us one of the four diagonal elements of the 
2-qubit density matrix, which are 

{\ + {q^zqz))±{{qiz) + {q2z)), {l~ {qizqz))±{{qiz) - {q2z)) 

Another operator which is particularly useful to us is 
the density operator. It is one of the most useful ways of 
going between the Schrodinger and Heisenberg pictures 
as in both cases we are dealing with the time-evolved 
operator, \'ip{t)){tp{t)\. Therefore, unlike ordinary oper- 
ators, there is no difference in the form of the density 
operator between the two pictures. The density operator 
is a proper vector in Hilbert-Schmidt space, although its 
evolution is different from other operators: 

p{t) = \m){m\ = c/i^(o))(^(o)ic/t = c/p(o)[/t 

The general form of a density operator in terms of the 
descriptors is 



P{t) 



where q'-"-' is the n-dimensional set of descriptors for the 
system, and {Pn} is the n-dimensional Pauli group. For 
a single qubit system. Pi = Ui (where i runs over four 
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indices and ctq = 1). For two qubits, Pi = ai® aj, and 
the qi{t) are therefore qij — qn ® q2j (again, qg — 1). 
That is, 

piit) = 

i 

Pi2{t) = ^{qii{t) q2j{t))ai (E) aj 



The two-qubit density operator can also be written 
p{t) = - f 1 (8) 1 + a.CT 1 + a (g) b.cr + ^ 



where 

ai = {qii{t)) ; h = ; *mn = {qim[t)q2n[t)) 

III. DIRECT CONSTRUCTION OF 
DEUTSCH-HAYDEN OPERATORS 

Using Deutsch and Hayden's method we can construct 
a density matrix from a set of qai gained by a certain 
preparation procedure. For example, consider a circuit 
of two qubits, with a Hadamard gate operating on the 
first and then a CNOT with qubit 1 as control and 2 as 
target. The qubits start at time f = in the zero state 

qi=(T0l q2 = l(8i(T 

The Hadamard gate takes them to 



qi = 



and then the CNOT leaves them as 





)1 \ 




/ 1$ 




-ay I 




q2 = j 












\ 1(5 





qi 




q2 




(1) 



The density matrix is therefore 

i 

which is the density matrix of the state |00) + 

A reasonable question to ask at this point would be: 
if we do not know (or do not care about) the prepara- 
tion procedure, can we then do the reverse operation, 
and construct a set of qai directly from a given density 
operator? 

Furthermore, how unique would such a construction 
be? In our example of |00) + |11), we can immediately 
see that the set ([1]) is not unique as we could swap the 
qubits in the preparation procedure. Performing H2 then 
CNOT2-^i would give us 



qi 




q2 



1 ®i ax 



(2) 



Under what circumstances in general do different prepa- 
ration procedures give rise to different sets of ^ai? Is 
this the only way that different sets of qai corresponding 
to the same Schrodinger state can be constructed? Es- 
sentially, how much can we trust what a given analysis 
of information flow tells us about the dependencies in a 
descriptor? 

We start by constructing a set of qai from a density 
matrix. The density matrix is dependent on the average 
values of the qai, so we must start by finding these. We 
can find the values of (qaiit)) by finding the components 
of p in the qi{0) directions, Tr{pqi). Such an operation 
can be simplified by noting that, writing p = J2n Pin ® 
P2n, we have 

Pv = Tr(^pi„ai (8)p2nO-j) 

n 

= ^'^r{pinCri)Tr{p2n(rj) 

n 

We therefore need only to look at the components in the 
sum where both traces are nonzero. The nonzero traces 
are: 

Tr(|0)(0|l)= Tr(|l)(l|l) =1 
Tr(|0)(0|a,)= -Tr(|l)(l|a,) =1 
Tr(|0)(l|a,)= Tr(|l)(0|(7,) =1 
Tr(|0)(l|fT,) = -Tr(|l)(0|a,) =i 



For example, consider the state |00) 
operator is 



1 11). The density 



P=|0}(0|®|0)(0| + |0}(1|®|0)(1| + |1)(0|®|1)(0| + |1)(1|®|1)(1| 



The nonzero components from this will be 1 ® 1, ax® ax, 
ay® ay and az®az, all of which are 1. The density matrix 



can therefore be written 



j=0 



which gives us 



{qu) = 
{q2^) = 
{qiiq2j) = % 

It is instructive now to look at the simplest possible set 
of qai that is consistent with these conditions (and also 
that qiy = qtxq-iz)- 



qi 




q2 




(3) 
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This, however, is not a well-formed set of qai as it does 
not form a basis in the Hilbert-Schmidt space of the two 
systems. A proper basis must have the following proper- 
ties (c/ I, pp3ff]): 

1. A''^ linearly independent operators: 3 {ay} such 
that ^ aij Cij — 

2. The operators of the basis are mutually orthogonal: 

3. The basis operators are complete: Tref^ = 1 

4. The operators must be Hermitian 

1. comes from the fact that the operators must span the 
entire space. 3 and 4 are consequences of the fact that 
all the operators are time-evolved from the t = Pauli 
operators, which are complete and Hermitian. This time 
evolution firstly cannot change the trace of an operator 
(3), and secondly is unitary so will take Hermitian oper- 
ators to Hermitian operators. 

We can see that the set ([3]) fails to meet the first of 
these criteria: it gives only 8 linearly independent oper- 
ators, rather than 4^ = 16. Contrast this with the qiiq2j 
operators obtained from ([1]), which comprise the entire 
set of Dirac operators, which are known to be linearly 
independent and space the entire space. 

If we now consider other qai which could give a well- 
formed basis, we can use these criteria to note also that 
qi and q2 cannot be identical: these co-ordinatise two 
subspaces of the whole space, and if the subspaces are 
isomorphic then the product of their bases will not give a 
basis for the whole space. Furthermore, no two qai can be 
identical, and neither can any two qiiq2j- In the first case, 
the resultant product operator would be 1 ® 1, which is 
already given by gio920 ilo = 1 is always understood), 
and hence in both instances there would be fewer than 
16 linearly independent operators. 

Let us see now how ^ may be altered to make it 
well-formed. The easiest way of making the combined 
operators span the system is to introduce ay instead of 
one (Tj. (preserving the average values) and for a 1: 



qi = I qixqiz 

a^,<Ty (g) l,a^ 



CTz, 1 (g) ay, ax 

q2 = I q2xq2z 

a^,ay® la^ 



Only two of these sets of these operators span the entire 
space: 



qi 



'z ^ 



-ay « ax 
ax'gl 




and 



1(g) ax 
qi = I ax (gay 



q2 



(Jz (» ax 

-ay g) ax I (4) 

ax g) t 



We now have three different sets of qai corresponding 
to the same Schrodinger state: (P), (0) and The 
interesting question now is to what extent such sets are 
unique - in the present case we have found three possible 
sets, are there any more that would represent the same 
state? In general, how many distinct sets would we be 
able to find? 

In order to answer this question we will first look at 
how many unknowns there are in the system (that is, how 
many unknown quantities are needed fully to determine 
a set of qai), and then consider how many constraints 
there are on their values. We start by writing a general 
form of a qi on the ath system: 



E 

k,l=0 



a^/cTfe g) ai 



(5) 



We know that qio = (720 = 1, so the general forms for 

qu g) q2j are {n, m G {x, y, z}) 

910920 = a ® a 
3 



91*920 = ^ Hi'^k ® ai 



3 

9io92i = E^™»^™®^" 

9ii92j = ^ a''i)kMj^lakajng>(Tia-„ 

We can see from this that if we can find the a^*;-* and bmn 
then we will have fully described the set of operators. 
aki for a given i and bmn for a given j each contain 16 
unknowns {k,l,m,n € {0,a;,y, z}). Now we know that 
9ay = qaxqaz, SO for cach systcm we have 2(16) = 32 
unknowns. Therefore, for our system here we will need 
64 numbers fully to determine a set of qai- 

So how many constraints do we have for the system? 
First, we have the conditions that no two qai are the 
same, and that they each have norm 1: 

Tr(9ii9ij) = 5ij 
Tr(92i92j) = 5ij 
Tr(9H92j) = 

where i,j G {Q,x,z}. Now, qao is fixed at 1 a, so in 
each of the above equations the case n = m = tells us 
nothing new. Therefore we have 3(9 — 1) = 24 constraint 
equations here. 

Next we have the constraint that each qai, i 7^ 0, must 
be traceless. To see why, consider two Hilbert-Schmidt 
spaces, one co-ordinatised by the Dirac basis {ai (g aj}, 
and the other by {91^ g)q2j}- There exists a linear trans- 
formation between the two, given by T. Because both 
spaces are of Hermitian operators, we require T to be 
unitary to preserve hermicity. Therefore, the elements of 
{9ii ^ 92j"} are a unitary transformation of the elements 
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of {(7i (8) CTj}. The trace operation is invariant under such 
transformation, so the elements qu q2j must have the 
same values for their trace as ai (Ei Uj - that is, 5^,5 j^. 
This gives us the four constraint equations 







Another set of constraints comes from the criterion 
that the elements of each are linearly independent. 
That is, that there exist some constants c„ for which 



E 



= 



(and for q2j)- Substituting ([5]) into this equation we have 

ikl 

We know that the Dirac operators are mutually linearly 
independent, so each (fc, I) term in the above sum must 
equal zero independently: 



aa 



"ikA 



As fc,Z G {0,x,?;,z} this gives us 16 equations for each 
system. However, we already know that a. 



0, 



so we have a total of 2(16) — 4 = 28 constraints here. 

It is interesting to note that up to this point, the con- 
straints on the values of qai have come exclusively from 
the structure of the Hilbert-Schmidt space - these are 
constraints on any physical system. We now move on 
to constrain our operators to a particular physical state. 
These constraints are the elements of the density matrix: 



Pi] = (9n'?2j) = Tr(gH(72j|0)(0|) 



(6) 



This gives us nine equations, but only eight constraints 
as we know (gio?2o} already (1). 

This gives us a total of 24 + 4 + 28 + 8 = 64 constraint 
equations for 64 unknown variables. Thus we see that, 
in general, the structure of the qai is fully determined by 
the physical system. That is, in general a given state 
defined by a density matrix has a unique repre- 
sentation in terms of Deutsch-Hayden operators. 

The caveat "in general" evidently applies as we have 
been dealing with a situation where the representation is 
not unique. Why is this the case for our above example 
00) + 1 11)? The reason is that we have a physical symme- 
try between the two systems, and this is picked up in the 
mathematics. All of the 'structural' constraints are, of 
necessity, symmetric between the systems and between 
basis states within each system (or else they would be 
making a physical statement). The physical symmetries 
are all contained within the structure of the density ma- 
trix. We can see from ^ that any symmetries of the 
density matrix will give rise to different sets of qai cor- 
responding to those symmetries. That is, if there is a 



symmetry of the density operator pertaining to the ele- 
ment pij^ then the same symmetry applied to qiiq2j will 
give us a physical state of affairs indistinguishable from 
the original. ^ tells us that such symmetries are the 
only means of generating separate sets of qaiS- 
Let us look at the density matrix in terms of the {ct^ ® 

3 

The symmetries that will give us different sets of qai are 
the symmetries of Oij : 

10 
10 
10 
1 

It is important to note at this point that we cannot 
have any symmetry that violates the constraints that 
we have imposed on the system above: so, for exam- 
ple, we cannot swap qix with another qai whose prod- 
uct with qiz is nonzero. Furthermore, qu cannot be 
swapped with qiiq2j unless j ~ 0: otherwise it would 
force q2j = 1 (g) 1 which would give less than 16 inde- 
pendent operators. To put that all in other words, we 
are restricted to the physical symmetry transformations 
of the system, non-physical transformations violating the 
constraints we have already laid on the system. In terms 
of symmetries of the density matrix, we are therefore re- 
stricted to swapping rows and columns. 
The symmetries of the density matrix are as follows ({} 
denotes a single transformation). Firstly, swapping single 
rows and the same numbered columns (ie reflections in 
the diagonal): 



{1 1 ^ 1 (g) 1} 



{ax ® 1 ^ 1 ( 



) a^} 



Then swapping two rows and two columns at a time: 



Ti ^ r2 and ci ^ C2 



1 (g) 1 ^ 1 0-3, 
10 1^0-2;®! 



n - rs and ci - C3 : { t ^ 1 ^ ay ® t 



Ti ta and ci <-> C4 



10 1^10 0-2 
10 1^0-201 



r2 - rs and C2 - C3 : { ^^^^^^^^l 



r2 ^ and C2 C4 
^3 ^ ^4 and C2 ^ C4 



1 CTj, ^ 1 (T^ 

CT2; 1 CTz 1 

1 CTy 1 

CTj, 1 ^ (72 1 



Any combinations of the above will also be symmetries 
of the matrix. 
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We can immediately discard symmetries of the form 
1 (g) 1 ^ 1 (g) : these will merely change the position 
of 1 1 in a set of qai, but we always have it as qao 
by convention. We can also discard several as qy is not 
independent of and qz - for example, swapping r2 
ra and C2 ^ C3 and swapping r2 ^ and C2 ^ C4 imply 
the same thing. We are therefore left with the following 
symmetries on a set of qai'- 



qiz 



q2x 



Qly ^ q2y 



qiz ^ 



qix 

q2x 



qix 

q2x 



liy 

I2y 



q2z 



liy 

<l2y 



qiz 
q2z 



(7) 
(8) 
(9) 

(10) 

(11) 

(12) 



We also have the combinations of the above (neglecting 
duplicates, and also leaving out the third in any given set 
of qi as it is implied by the first two): 



qix 

qiz 



^ q2x 
^ q2z 



q2x 
qix 



qiz 

q2z 



q2x 
qix 



q2z 
qiz 



^ qiy 

^ q2y 



qiy 
q2y 



(13) 



(14) 



(15) 



(16) 



These are all the symmetries present: if we take further 
combinations of (|13m6[) with ([7l fT2|) we get a symmetry 
already given. Therefore, the identity transformation 
qai qai plus the transformations ([7|) - p6|) form the 
group of relevant symmetries of the density matrix as 
applied to the qai operators. 

We can now use this group to generate the set of sets 
of qai that will correspond to the state |00) + |11). Now, 
because the only way to generate a different set of qai is 
by one of these symmetry transformations, each element 
in the set of sets of qai will have orbit 1 when acted on by 
the group of symmetries (we denote this by Gp) - that 
is, any element can be reached by any other using one of 



the group elements. For example, ([2]) is generated from 
(P) by applying dTO]), and (g]) by applying We can 
therefore use ([T]) to generate the entire set of allowable 
sets of qai'- 



(2) 



(3) 



(4) 

qi = 



(5) 

qi = 



(6) 

qi 



(7) 

qi 



(8) 

qi 



(9) 

qi 



(10) 

qi = 



(11) 

qi = 



a, 

—a. 



1 ® 

(7x (Jx 

CTj; ® 1 



1® <7y 

CTx ® CFy 



<7x 

CTx ® CFy 

'y 



CTx 

a,, 



y 



-a-y a-x 
az ® (Jx 

CTj; ® 1 



az ® CFx 

CTj; ® 1 

-ay ® ax 



1 (g) tJa; 

CTj, ® ay 





(2) 

q2 



(3) 

q2 



(5) 

q2 



q^^ 



(8) 

q2 



qf 



(11) 

q2 



1 ® (Tj^ 




—a,, 



Cx ® CTz 








(Ta; eg) 1 



51 
' Cfx 



Or, 



-Oy 

az ® a^ 
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1 ® CTx 

Qi =1 cr^(S)a^ 



(12) 

q2 




There is one final question that needs to be answered 

about these alternative sets of qaV- do they give the cor- 
rect results when they subject to further evolution of the 
system? Do all these sets give the same average values 
after arbitrary unitary evolution? 

It is fairly straightforward to prove that this is the 
case. Consider two set of operators sharing the same 
average values, qai and q'^^. qai has been constructed from 
a circuit and so is known to give the correct evolution, 
ie p{ti) = Y.ij{<lu{'^)<l2]{f))(7t ^ cFj. Because qai and g^ 
share the same average values, we can write 



ij 

Let U be the transformation taking place between t and 
pih) = Up{t)U^ 

ij ab 

We showed above how any basis can be written as TV„ (8) 
Pit^) = J2(U''{t-to)ai^ajU{t-to)) 



ab 



Now 



Tr(cro (8) CTf, Uai ® ctjU^) = Tr(C/[/Vo ®abUai® (JjU^) 

= Tr(?7Vo (g) (Tb Uai <8> (Tj) 

_ j 1 if (Tj (8) (Tj = U^aa <Sl (Th 



otherwise 
Therefore 

(Ta 

ab 

= X!(?7t(<i - to)(7a (8 (JbU{ti - to)) a a ® ab 

ab 

= - t)q[a{t)q'2b{t)U{h - t)) aa ® ab 

ab 

= ^{l'laiil)(]'2b{tl))'^a®Ob 
ah 

In other words, for arbitrary evolution, the evolved den- 
sity operator can be written in terms of the evolved q' 

operators, and 

(au(ii)<?2,(ii)> = {<iu{ti)<l2j{ti)) 



We have therefore shown that in the Dcutsch-Hayden 
picture a given Schrodingcr state is represented by the 
Gp-set of sets of qai which is acted on by the group Gp 
of symmetries of the density matrix p of the state. Each 
element of the Gp-set has orbit 1, and for all practical 
purposes any element may be used in place of any other. 
We therefore see that in an analysis of information flow, 
the possible descriptors for the systems under considera- 
tion differ only by the physical symmetries of the systems 
- that is, the description of information flow and depen- 
dencies will not differ in physically significant ways be- 
tween different sets of descriptors for the same situation. 
Any set that we chose will give us the correct description. 



IV. THE MEASUREMENT INTERACTION 

We now turn to measurement-type interactions in 
terms of the Deutsch-Hayden picture. The simplest form 
of a measurement interaction is full decoherence, where 
the system is measured by an ancillary system of the 
same size to which we do not have access. The simplest 
model of this interaction is that there is a CNOT gate 
with the system as the control and the ancilla as target. 
For example, consider a one-qubit system and a one-qubit 
ancilla. Without loss of generality, we can take the state 
of the ancilla before interaction to be |0} Qa = 1 (8 cr. 
The state of the system before is given by (\s{t). The 
state of the system afterwards corresponds to 

qsx{t) (8 a^ 
= I qsv{t) (8 ax 
qsz{t)^l 

That is, the only nonzero (qsi) is qsz- Now, the diagonal 
elements of the density matrix are 1 ± cr^ - so what we 
have is the standard action of decoherence, where the off- 
diagonal elements of the density matrix written in the 
decoherence basis go to zero and the diagonal elements 
are unchanged. 
U We can extend this to any size of system if we have an 
ancilla at least as large as the system being decohered. 
Then each qubit of the system will perform a CNOT 
with a given qubit of the ancilla, leaving only the {qnz) 
as nonzero. In any size system, the diagonal elements are 
given by the coefficients of combinations of 1 and CTz 's - 
that is, combinations of (1) and (q^) only. Thus again 
we have the off-diagonal elements becoming zero and the 
diagonal elements unchanged. 

Moving away from decoherence, let us now suppose 
that we have access to the ancilla. In our two-system 
example, the operators for this ancilla will be 

1 (8) (Ta; 

qa = I qszit)'S>ay 
qszit) (8(T2 

We can see from this that only the information contained 
within the qsz{t) operator is picked up by the ancilla. 
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This is because, were we at a later time to have access 
only to the ancilla, the information about the first system 
that could be extracted is qsz{t) only. Again, this is what 
we want: the qsz{i) components give the probabilities of 
outcomes of measurement in the measurement basis (as 
they give the diagonal elements of the density matrix). 
This is all the information about a state that can be 
gained by measuring it. 

We have seen what happens when we perform a mea- 
surement in the |0),|1) basis. What if we wish to measure 
in a different basis? In such a case, the system would 
still only have nonzero {qz) component (and the ancilla 
only pick up information from qz) - but after the qai 
had been rotated in the change of basis. For example, 
suppose we wished to measure the system in the basis 
0) ± To change to this basis from |0), |1) we perform 
a Hadamard transformation. The system operators are 
therefore {az,—<JyT<Jx)- Performing a CNOT with the 
ancilla then gives us 

(qszit) (E) CTx \ 
-qsy{t) (Tx 
qsx{t)(E)t J 

That is, the coefficient of the non-zero component cfz is 

(qsxit)). 

What we have in effect done here is change the meaning 
of a as the basis changes, and compensate by changing 
the coefficients, a becomes with respect to the new basis 
rather than |0), |1) and the coefficients will necessarily 
change. 

In modelling the measurement of the system as a 
CNOT gate with an ancilla as the target, we have been 
dealing solely with projective measurements. For com- 
pleteness, we would like to look at the case of generalised 
measurement, where the states of the measurement basis 
are not necessarily mutually orthogonal. 

A POVM measurement can always be modelled by 
adding a second system to the one being measured, per- 
forming a projective measurement on both systems and 
then tracing out the second system. Now, we have seen 
that the density operator in any basis can always be writ- 
ten as 

l(g)U-ha.(7(8)U-hll(8)b.cr-H^tm„crm(g)cr„ (17) 

mn 

When the second system is traced out this becomes 
a + a.a 

All the interesting information is therefore contained in 
a = qiz- 

Take the single-qubit system we measured above. We 
add a second system, choosing it to be in the |0) state. 
Now we measure the joint system in the Bell basis. In 
terms of the qai this means performing a Hadamard then 
a CNOT gate. At the end of this procedure the operators 



for the two systems are: 

/ qiz{t)(^a^ \ / l®(Tx \ 

qi = -qiy{t)(E)(Jy q2 = qiy{t)®(Jy 

\ qix®^ j \ qixit) ®az j 

To complete the measurement, we must add a two-qubit 
ancilla and perform CNOTs between the two systems and 
the ancilla. We saw above that this will make all the av- 
erages for the two qubits individually zero except {qz)- 
Therefore, in the generalised as in the projective mea- 
surement, what can be gained at measurement is knowl- 
edge of qsz ■ 

In what way, then, is the generalised measurement 
more powerful than a projective measurement? In terms 
of the Deutsch-Hayden operators, the difference comes 
in the choice of the basis for measurement. There are a 
much greater range of operations that can be carried out 
on two systems rather than one, and so there is more op- 
portunity to as it were 'move' different information into 
the q^ 'slot' from where it can be picked up by the mea- 
surement. 

V. RELATIVE STATES 

We have seen what happens to the state of a system 
which is measured, and where the ancilla which per- 
formed the measurement is subsequently traced out. We 
now look at what happens when the ancilla is not ignored, 
but rather one state of it is singled out - in such a case, 
what is the relative state of the system that was mea- 
sured, as given by the Deutsch-Hayden descriptors? The 
above results will, of course, be the special case where 
the relative state of the ancilla is completely unknown - 
that is, 1. 

As is well known, if we have two systems with joint 
density matrix pi2 and a state of the second system \(3){l3\ 
(this can be one of the elements of a POVM), then the 
state of the first system is given by the partial trace 

Pi=Tri(pi2|/3)(/3|) 

Writing the density matrix as (|17p we have 

pi = l.(/3|]l|/3)+a.a.(/3|l|/?) + IL.(/3|b.a|/3) 

+ y^JmnCr7n-{P\crn\P) 

mn 

= l(l + (^|b.a|^))+a.a 

+ y^Jmn{P\0'n\l3)0'm 
mn 

We want to write this in the form o-'iO'i, which gives 
us 

= + y^Ji„ {l3\an\P) 

n 

In terms of the qai we therefore have 

(qii) — > {qii) + C^qiiq2n{f3\crn\f3)) 

n 
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This can be achieved if the q[ for the first system relative 
to the state |/?)(/3| of the second system are 



This expression is useful if we have the relative 
Schrodinger state. We can also frame this is terms of 
relative Deutsch-Hayden operators. If are the opera- 
tors for the second system when it is in the state \(3) then 
we know that 

(0-«)|/3)(/3| = {q'2n) 

(ie the coefficient of cr„ in the expansion of the density 
operator). The operators for system 1, q[^, relative to 
the operators (/jj of the second system, are (unprimed 
operators are for the original state of the systems) 



9ii = 9iHl + y^('?2n> 12, 



(18) 



We check this for the case where the second system is 
simply discarded, and its density matrix is therefore 1. 
In that case we have (remembering that n runs from 1 
to 3) 



q'li 



= qii 



1 



q2n 



= qii 

which is as we found above. 

We can now use the expression for relative states to 
look in further detail at the measurement interaction. Let 
us consider here a specific example, where the system to 
be measured is in the state |0) + |1), with descriptors 



qi 



a 



1 



Now we add the ancilla and perform the CNOT opera- 
tion: 



qi = 




q2 = 




(19) 



Now we can look at the states of the first system relative 
to those of the ancilla. We choose the relative states of 
the ancilla to be in the computational basis, and use (fT5)) 
to give us the relative state of the system. The relative 
operators of the system are therefore 



qi 



—a,. 



(20) 



q'l' 



) a 



(21) 



a 



where the primed are relative to |0)(0| and the double- 
primed to |1)(1|. We can see from these qai that the 
states of the systems in both cases are eigenstates of one 
of the z± operators (recall that these operators give the 
probabilities for and 1 in the computational basis) . An 
interesting question at this point is to ask what the cor- 
responding operators for the ancilla are - that is, what 
are the complete sets of primed and double-primed oper- 
ators? 

To answer this question we must pay close attention 
to the physical situation of the measurement. We have 
in a rather cavalier fashion introduced the states |0)(0| 
and |1)(1| of the ancilla, which is part of the joint sys- 
tem. Where do these states come from? What we are 
always meaning in a situation such as this is that the 
system under consideration (here, the ancilla) has been 
measured in a basis of which the relevant state is a part, 
and then one basis state over others 'picked out'. Physi- 
cally, it is picked out by being the state which is relative 
to some other state that we are interested in. At some 
point, then, in a tractable analysis we will have to simply 
stop and pick out the state relative to which we wish to 
find other states by fiat. At this point we construct the 
operators for this 'ultimate' state from the density oper- 
ator that we wish it to have, and then find everything 
else relative to it. Where this operation takes place is 
entirely a matter of convenience. What must be remem- 
bered, though, is that at this cut-off point we then lose 
the record of the interactions of that system with other 
systems. For example, we could not in the current sit- 
uation choose to construct the operators for the ancilla 
from the ground up - that would give 



1. 



q2 = I ±1 
±1 



comparing this with p9)) . we see that we have lost the 
information about the first system that was contained 
in the ancilla. 

In this situation, what we need to do is push the ulti- 
mate state back one system further, and have the ancilla 
measured by a third system. If we start that system in 
the |0) and perform CNOT with the ancilla as control 
and the third system as target then we have 



q2 




qs 




The operators of the ancilla relative to the states |0) and 
|1) of the third system are, respectively, and q^ : 



q^ 



q2(l ±93^) 
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((1 (g) CTa; ± (Tx ® (Jz) ® {(Jx ± \ 
{(Jx®(Jy±CTx®(Tz)®{(7x±(yz) (22) 
{ax ®az±crx® <Tz) ® (ctj: ± c^) / 

These then are the operators for the anciUa to be in the 
states |0)(0| and |1)(1|, relative to which the first system 
states are ^ and ([2T|). 

If we look at , (|2ip and (I22p we notice something 
interesting: the sum of the relative qai is (neglecting a 
factor of 2) the original qai for the system. This is, in fact, 
the case for any size system where the states relative to 
it are a complete POVM ~ the sum of the qai relative to 
each of the states of the POVM is the qai for the system 
before measurement. The proof is straightforward. Using 
(fTS]) we have 

m 

i=l n i 

Now 

Y,{a\A\a) = Y.MA\a){a\) 

a a 

If {\a}{a\} forms a POVM then J2a = 1> so 

J2{a\A\a) = JrA 

a 

Therefore, as the an have zero trace, we have 
^q^i^ = qi(m]l + ^Tr(cr„)92«) 

i=l n 

= mqi 

Physically, then, we have performed a POVM mea- 
surement on the ancilla system, and are now looking at 
states of the first system with respect to elements of the 
POVM. 



VI. MIXED STATES 

So far we have been considering only the representa- 
tion of pure states in the Deutsch-Hayden picture. We 
shall see that one of the advantages of this picture is 
that the representation of mixed states does not differ in 
kind from the pure-state representations. We will also see 
that the Deutsch-Hayden representation of mixed states 
is transparent to their physical origins. 

The Deutsch-Hayden operators qai form a basis in the 
Hilbert-Schmidt space. This is not, however, an arbitrary 
basis: the elements of the individual combine to form 
the elements of the basis. We have, in fact, a product 
basis on the Hilbert space. We can therefore consider 



the space of the system as a product Hilbert space, of 
factor spaces both of which (or each of which for more 
than two systems) are co-ordinatised by Qq . These factor 
spaces will, of course, evolve over time - the product does 
not remain fixed. However, at each time t the space of 
the system may be written as the product Hi {t) ® H2 (t) 
where {qii{t)} forms a basis on Ti.i and {q2j{t)} forms 
one on Ti.2- 

In order to see to what these factor spaces corre- 
spond, consider the situation at time t = 0. In this case 
qi = cr (g) 1 and q2 = 1 (g) cr, and in the Schrodinger 
representation we have the state |00). The spaces co- 
ordinatised by the qa are therefore transparently the 
spaces corresponding to operators of the two systems. 
Let us look at, for example, an operator that operates 
only on the first system: (g) 1. This can be written 

1 (g) 1 = ^ a^CTi ® a 

i 

= atqiijO) 

i 

Therefore at a given subsequent time t we have 
A{t) ^^a,qii{t) 

i 

That is, anything that can be said about the first system 
alone is contained within the qu alone, with no reference 
to other qai- The space co ordinatised by quit) at any 
given time is therefore the space of operators on the first 
system alone. Thus, at any given time, the overall space 
of the system is a product space of the individual system 
spaces. Unlike in a standard Hilbert space representa- 
tion, this is the case always, regardless of whether the 
systems are entangled or not. 

In none of the foregoing have we had to say whether 
the qai correspond to pure or to mixed states. It might be 
questioned whether mixed states can indeed be accom- 
modated in this picture - after all, we have constructed 
qai by evolving them unitarily from a fixed state, and we 
know that unitary evolution alone cannot produce mixed 
states from pure. 

The key to representing mixed states is the fact that 
all mixed states can be considered as pure states where 
one or more systems have been traced out [1, ppllOff]. 
The Deutsch-Hayden operators for a mixed state on a 
system a will be the qai which have been evolved as part 
of a larger, pure system. For example, we have found the 
qai corresponding to |00) + |11) on two systems. If we 
trace out, for example, the second system then we are 
left with the maximally mixed state |0)(0| + |1)(1|. This 
is represented in Deutsch-Hayden terms by 

(<Jz(E)ax \ 
-ay (g) ax (23) 
axfEit J 

We can see that one of the most immediate differences 
between pure and mixed states is the dimension of the 
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non-trivial part of the subspace for the system. A pure 
state of the first system would be represented by opera- 
tors of the form (g) 1 - for example, |0) -I- |1) is given 
by 



qi 



by contrast, the subspace of a mixed-state system is non- 
trivial in more dimensions, corresponding to the systems 
which have been traced over to get from the pure to the 
mixed state. Note that in the Deutsch-Hayden picture 
there is no such thing as tracing out a system: whether 
or not we are looking at the second system makes no dif- 
ference to the way in which operators on the first system 
are constructed. 

It is a necessary and sufficient condition of a state being 
mixed that its qai cannot be written non-trivially on the 
same size subspace that a pure state of the same system 
would require. If this were possible, then the operators 
in the non-trivial part of the subspace could be written 



= [/V ® a ( 



>crC/ 



for some unitary operator \J . The average values of the 
q could then be written 



(q) = Tr(q|0)(0|) 

= Tr((T ®a®. 
= Tr((T ®a®. 



.®aV |0)(0|) 
10- C/|0)(0|t/1') 

1 cr p) 



Note that p in the final line is pure, which gives the con- 
tradiction: these are the average values of a basis set in 
some pure state, which cannot give the same outcomes 
as a mixed state. We will therefore always need a larger 
subspace to describe a mixed system than a pure one. 

As well as distinguishing pure and mixed states in this 
way, we can also frame the standard condition Trp^ < 1 
in terms of the qai (here for 2 systems): 



Pl2 



1 

16 



("ooTr(l «) 1) -t- ^ aojao„Tr(l ® CTjCr„ 

a|,Tr(a2 



ijmn 



+ aioamoTr{aiam 1) 



)) 



1(4 + 4(V«' 



ij 



1 ^ 



The condition for a mixed state is therefore 

3 

E(<ZH)' + (92,>' + (gi.g2,)' <3 (24) 

The Schmidt decomposition condition for pure states 
can also be framed in terms of the Qai- The Schmidt 
decomposition of a pure state density matrix is 



PAB = Yaij\i){j\ ® 1001 



Writing this in terms of cr-matrices, that is 

PAB = a{l + (Jz) ® (1 + (Tz) + h{ax - lay) ® {a^ - icTy) 

+ c{ax + icTy) (ctx + io-y) + d{t ~ az) §5 (1 - (Tz) 

The coefficients are therefore (in terms of the qai) 



(10 1) + {qiz) + {q2z) + {qizq2z) 

{qixq2x) - {iqiyq2x) - {iqixq2y) ~ {qiyq2y) 
{qixq2x) + (iqiyq2x) + {iqixq2y) - {qiyq2y) 

(1 1) - (qiz) - {q2z) + {qizq2z) 



These follow the usual Schmidt rule + 



(f = 1, 



which when expanded in these terms gives (|24p . 

As well as distinguishing mixed states when we are 
presented with them, we can also construct the qai for 
mixed states in the same way that we constructed them 
for pure states. We can either construct the purified state 
using a circuit, or we can directly construct from the 
density operator. Let us look at both of these in turn for 
the simple states |0)(0| + |1)(1| given by (^5)) . 

Firstly we will try the purification method. This is 
straightforward [1, ppllOff]: take the mixed density op- 
erator in its diagonal basis and then add another system 
with the same basis, to make a pure state whose Schmidt 
decomposition for the original system gives the basis in 
which the density operator is diagonal. In the case of 
|0)(0| + |1)(1| we add a second system with basis {0, 1}. 
Possible purifications are |00} + and |01) -1- |10). We 
now find the corresponding qai for the whole system, and 
then take only the qai for the original system. In the first 
case we would have ([^5)). The second can be found by 
flipping the second bit in a state |00) -t- |11), which has 
the effect on the qai of being acted on by 't ® ax- As it 
happens, this gives us back the same qai^ ([23]) . 

Let us now try direct construction. We have the den- 
sity operator for the state. 



P 



1 
1 



which tells us 



ii=l 



{q^) = 



12 



We know that we will not be able to express the qi as 
combinations of a (in this case as well this is obvious), 
so we will add an extra two dimensions to the subspace 
under consideration and try combinations of ct^ ® cjj . We 
have quite a lot more freedom in this construction than 
in that of the pure states, as we do not care what the 
operators for the ancillary system are. All we need to 
make sure is that the operators for the system we do 
care about co-ordinatise a well-formed subspace of space 
for the whole system. It is easiest to start off with a in 
the first position, and then add where needed, giving 

qi = CTy^CTx 

We verify that this is indeed a representation of |0)(0| + 

|1)(1| as they are the qi operators from the set q^^ in 
the previous section. 

This representation of mixed states is a very physical 
one. The Deutsch-Hayden picture does not recognise any 
form of evolution other than unitary evolution, and this 
is reflected in the fact that even when a set of Qai for 
a mixed state are directly constructed as above, their 
structure shows that they are in fact part of a larger 
pure state. The size of the mixed-state operators can be 
found from the purification method used above. An extra 
system is needed to purify each system that is mixed ~ 
so the Deutsch-Hayden operators will in general be twice 
as long for mixed states than for pure states of the same 
sized system. That is, a mixed state of N qubits will in 
general be part of a pure state of dimension 2.2^, and 
the Deutsch-Hayden descriptors will have length 2N. 

VII. REDUCED DESCRIPTORS 

We have seen that an n-dimensional mixed system 
needs to be written on a larger space than a correspond- 
ing pure system. We will now look in detail at the ability 
to write a descriptor on a reduced space, in particular 
with reference to entanglement between systems. 

Let us take a three qubit system, the overall state of 
which is pure, and then neglect the third qubit (that 
is, we are only interested in the descriptors for the first 
two systems). If the first two systems together are not 
entangled with the third then their descriptors will satisfy 

Pi23 — P12P3 
(qi23)a' ® ® = (qi2)a' ® cr (q3)cr 

(qiq2q3> = (qiq2)(q3) 

Writing this out in full, with U giving the evolution of 
all three systems gives us 

Tr(|0)(0|J7tcr(8)cr(8)crt/) 

= Tr(|0)(0|i7V (g) (g) lU)Tr{\0){0\UH (g) Kg) aU) 
Tr((7 (g) cr «) CT U''\0){0\U) 

= Tr(|0)(0|C/V(gcr(glC/)Tr(l(g)l®cr U^O){0\U) 



Writing a ® b (g) c = W\0){0\U we have 

Tr(cr g)crg)cr ag)bg)c) 

= Tr(|0) (0| C/V g) (T g) lU) Tr(l g)l(gcrag)b(gc) 

= Tr(|0)(0| ijt tJ g) cr g) lU) Tr(a (g b g) ac) 
Tr(o- (gCT(gla(gbg)Il) Tr((Tc) 

= Tr(|0)(0| J/V g) cr (g IC/) Tr(a g) b g) l)Tr{ac) 
Tr(o- (g(T(gla(gbg)IL) 

= Tr(|0)(0| [/V (g cr (g IC/) Tr(a ® b g) 1) 

Now in general a, b and c are sums, X^i^o^^"^^' "^ti^re 
t'^'^ t'^'' t^Q^' is the 1 g) 1 component of a (g b g) c. As this 

component is never zero, ip"''''^'* 7^ 0. Therefore the trace 
of a, b and c will always be 1 (neglecting normalisation) 
- and hence Tr(ag)bg) 1) is a simple, constant numerical 
factor which can be neglected here, and can be incorpo- 
rated into the normalisation of the descriptors in the end. 
We therefore have 

Tr(|0)(0| J/V (g cr (g IC/) = Tr(cr g) cr g) a a g) b g) 1) 
(qiq2) = Tr(cr g) cr a g) b) 

Now let us define U12 as the elements of the matrix U 
that act only on 7ii2. We then have 

(qiq2) = Mag)a Ul^\0){0\Ui2) 

(qiq2) = Tr(C/iV®(7[/i2 |0)(0|) (25) 

What then are the expressions Ul20' g) cr C/12 exactly? 
They are the same as the descriptors C/V(gcr(gl?7 where 
only the action on W12 is taken into account. We can 
therefore see that they are the "simply reduced" qiq2; 
that is, the descriptors with only their H12 components. 
For example, if qi = Ox®o z®(Jy then the simply reduced 
form on Jin is o^® (^z- 

At this point we need to introduce some terminology. We 
will say that: 

A descriptor q is said to be represented 
by an expression n when it is the case that 
(q) = (n) 

We saw examples of representation previously, where 
many different sets of q's described the same density ma- 
trix; each set was represented by the others. However, 
representation is not restricted to such examples on the 
same Hilbert space: in general the dimensionality of n 
and q differ. 

Looking at (|25p we can see that this is a case of rep- 
resentation. We have shown that qiq2 S ^123 can be 
represented by vX^a (g crUn G H-n- This is therefore the 
key to representing pure states. When looking at a sys- 
tem s, only the elements of the descriptor on Tig need be 
used. When considering only that system, the compo- 
nents on other spaces make no difference to the average 
values, and might as well be replaced by 1. 
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As a consequence, when looking at the past evolution 
of the system we may as well replace the action of U on 
H3 by 1 as well. 

(It becomes necessary here to introduce some notation. 
There is the potential for confusion between descriptors 
with different numbers of components. When the context 
does not make clear on which space they have support 
then we will write for example [qi]i2 for the descriptor 
qi with components only on H12 - that is, Ul2(J®t-Ui2-) 

Another way of seeing this is by considering the fol- 
lowing way of expressing the descriptors for the first two 
systems: 

(qiq2) = Tr(|0)(0| L/V0cr0lC/) 
= Tr{U^\0)(0\U a(S)a(E)l) 
= Tr(J7|2|00)(00|t/i2 cr® cr)Tr([/]|0)(0|J7) 
= ([qiq2]i2) Tr(C/]|0)(0|C/3) (26) 

So what happens when there is no entanglement with 
qubit 3 is that the expression Tr([/3 |0)(0|f73) becomes 
irrelevant; it is not zero when ([qiq2]i2) is not zero. 

We have seen what happens when there is no entan- 
glement between the two systems. What changes then 
if the first two qubits are in a mixed state - if they are 
entangled with the third qubit? If we restrict ourselves 
to simply that system, what can we say about the qubits 
and their descriptors? The most obvious point is that we 
will no longer be able to describe them on 'H12 using the 
simply-reduced descriptors. The information contained 
in 7^3 becomes relevant to describing the past evolution 
of the two-qubit system - as shown in Chapter 2, we 
require 7ii23 to write them on. 

We now have a problem. As can be seen from all the 
foregoing work on the Deutsch-Hayden approach there 
is no operation in the formalism that corresponds to the 
action of tracing out a system, which is the usual way 
of getting a mixed state representation. Subsystems and 
their descriptors have been considered in their physical 
context, as part of the wider system, and their descriptors 
written on the space of that wider system. The question 
is: given that there is no physical operation correspond- 
ing to tracing out a system, what happens when we have 
no knowledge of anything other than the system under 
consideration? The solution is prompted by the foregoing 
paragraph: that ignoring the action of the third system 
means ignoring Ti.^. So if we want to restrict our atten- 
tion simply to the system of qubits 1 and 2 then we can 
only have access to 7ii2 on which to write our descriptors. 
Hence the problem. 

This may appear to be labouring the point somewhat 
as it is obvious what the solution to this is going to be. It 
is, however, instructive to show exactly how that comes 
about in the Deutsch-Hayden approach, and why it is 
needed - why we cannot just deal with the descriptors as 
we have them on 7^123. 

The solution is that we have found a set of operators 
on 7^123 which we know can be simply reduced to ?^i2. 



If we can write our 'mixed' descriptors in terms of these 
'pure' descriptors then we have a way of representing the 
descriptors associated with mixed states on 7^12. There 
are two ways in which we could do this. Either a com- 
bination of the descriptors associated with pure states 
could represent the descriptor we are concerned with, or 
else the pure descriptors could form a complete basis on 
7^123 for all possible uV (g) cr (g) lU. 

We know from standard quantum mechanics that the 
first is going to be the case - this is the standard way of 
writing mixed states as a mixture of pure ones. We can 
show this in the present case. Let us look again at (|26p . 
We can re- write this as 

(qiq2) - Tr(|00>(00| {U^a aUi2Mul\0) (OlU^))) 

Writing 

n m = Ul^a ^ CT[/i2Tr([/]|0)(0|i73) 

we have 

(qiq2) =Tr(|00)(00| (n®m)) 

That is, we can represent qiq2 on 7^i2, by the expression 
n g) m. Now, n (g) m can be expressed in terms of a 
complete basis {Ai} € 7^12: 

(qiq2)=Tr(|00)(00| (J^A^A.)) 

i 

We know one candidate for {A^}: the set of all possible 
descriptors on 7^i2, namely C/V (g) all. We can therefore 
write 

(qiq2) = Tr|^|0)(0| w.Uj a ® aU^j 

= Tr^|0)(0| J2w,Ulc7®c7®lu}j (27) 
= ^Wi{c[uq2i) (28) 

i 

where 

= MUja ® aU^ Ul^a ® aUu) Tr(C/]|0)(0|[/3) 

and where the q„i (given by the Ui) are those descriptors 
on 7^123 that describe qubits 1 and 2, unentangled with 
qubit 3. 

What about the second possibility, that any given 
qiq2 G 7ii23 can be written as a mixture of the pure 
qiiq2i, rather than simply being represented by them? 
This would require us to be able to write 

qiq2 = ^A,C//cr(gcr® If/, (29) 

i 

where again the Ui give qubits 1 and 2 unentangled with 
qubit 3. We have written A^ here rather than Wi as they 
do not have the same form: 

A, = Tr(C//(T g) (T g) IfJ, Ul^^a g) ct g) i;7i23) 
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where U123 gives the first two qubits in a mixed state. 

We can show if this is the correct form of qiq2 by 
considering the physical predictions that it would give. 
We know that the correct values for the averages of these 
descriptors are given by ([28)1 . Writing qiq2 as in ([29]) 
would give as the average values 



(qiq2) -Tr |0)(0| 



Comparing this with ((27)) we see that these can only be 
identical if Xi = Wi. Suppose that this is the case. We 
would then have 

Mula ® aU^ Ul^a crt/12) Tr(/7]|0)(0|C/3) 
= Tr{Uja (S)al»lU^ Ul^^a (g) cr (g) aC/123) 

Now let us write 

U^a (g) CT ® aC/i = Uja g) alii g) 7^ 

[//jgO- (g) cr (g) i;7i23 = C^iV ® cr^l2 ^ C 

This gives us 

MU^^a^aU^ Ul^a^aUi2) Tr(C/j|0)(0|C/3) 
= Tr([//cr g) crC/, Ul^a (g crL/12) Tr(7ic) 

which would require that for all i such that Tr{uja g) 
aU., Ul^a®uUi2) + 0, 

Tr([/t|0}(0|{/3)=Tr(7,c) 

which is a contradiction, as the left hand side is a con- 
stant value for all i whereas the right hand side is a vari- 
able. 

It is therefore not the case that qiq2 itself can be writ- 
ten as a mixture of descriptors for pure states, but it is 
the case that it can be represented by such. 

This is a very interesting situation. What this shows 
is that in the Deutsch-Hayden approach the term "mix- 
ture" is fundamentally incorrect. The descriptors for 
such systems can be represented by a mixture of descrip- 
tors for pure states, but they are not themselves identical 
with a mixture, even when looking at the pure descrip- 
tors spanning the original Hilbert space. All of this is 
very different from a standard density-matrix formalism, 
where the pure state density matrices co-ordinatise the 
space of possible density matrices and are hence 'special' 
in a way that the mixed density matrices are not. In- 
deed, the mixed density matrices are given as different 
types of entities from the pure state ones (a fact empha- 
sised by the standard Schrodinger notation, which cannot 
deal with them). 

In terms of the Deutsch-Hayden approach, mixed and 
pure states are both described by sets of descriptors on 
a Hilbert space of certain dimensions, this corresponding 
to the number of systems involved in the evolution of the 
system under consideration. The difference between pure 



and mixed comes when we wish to reduce the Hilbert 
space on which we write our descriptors - that is, when 
we wish to ignore certain subsystems in the analogous 
operation to tracing out systems. Only then does any 
notational difference arise, some systems needing a sum 
of descriptors to describe them, some able to use the 
reduced descriptors. However, they are still both on a par 
with each other: the descriptors corresponding to pure 
states are not 'special' in any fundamental way, and they 
do not co-ordinatise the space of all descriptors. Each 
are as fundamental as the other, without the definition 
of one being dependent on the other. 



VIII. ENTANGLEMENT SWAPPING 

We will now look at the flow information in an entan- 
glement swapping situation. In the Schrodinger repre- 
sentation, the particular example that wc will use is the 
following Q. We have four qubits in the overall state 
1$): 

1$) = (|00) + |ll))i2®(|00) + |ll))34 

= |00)23lOO)i4 + |01)23|01)i4 + |10}23|10)l4 + |ll)23|ll)l4 
= (|00) + |11))23 ® (|00) + |ll))i4 

+ (|00)-|11))23 g) (|00)-|ll))i4 
+ (|01) + |10))23 g) (|01) + |10))i4 
+ (|01)-|10»23 g) (|01)-|10))i4 

If we then measure qubits 2 and 3 in the Bell basis, we 
can see that qubits 1 and 4 become entangled, and the 
measurement destroys the entanglement within the pairs 
(1,2) and (3,4) - it has "swapped" to the pair (1,4). 

The locality issues of this operation are well known: 
the qubit pairs (1,4) and (2,3) could be sent to opposite 
ends of the galaxy after the original entangling opera- 
tions, and yet a subsequent operation on (2,3) will entan- 
gle (1,4). Even worse, qubits 1 and 4 could be separated 
before the measurement on (2,3), yet after it they are in 
a maximally entangled state, which can for example be 
used to teleport between them. 

We start our analysis with the four qubits in the state 
1$): 



qi = 



qa 




q2 



q4 




We now rotate qubits 2 and 3 to the Bell basis by per- 
forming a Bell gate, giving 



qi 




q2 
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1 (X) 1 (g) CTa; (g) 1 

q3 — I —1 €5 (Ta: (8) (Ty (g) 

1 (g) CTa; (g) (Tz (g) (7:1 



q4 





Now we perform the measurement on the pair (2,3) in 
the Bell basis by introducing two further qubits and per- 
forming two CNOT operations between (3,5) and (2,6), 
leaving us with the complete set of descriptors 



qi 



q4 



qs 



qe 



q2 



qs 



If we now look at these descriptors, we can see some 
very interesting things. First, let us look at the density 
matrices for the pairs (2,3) and (1,4). For both pairs, 
Pab = 1 g) 1 and pa,b — 1, from which we can see that 
the pairs are not entangled. Next, we note that neither 
pair is in a pure state, as we cannot reduce their de- 
scriptors simply to the elements on the Hilbert spaces 
for those systems. Furthermore, if we look at the pairs 
(1,2) and (3,4) which were originally entangled, we see 
that they too have pab = Ig)!, Pa,b = 1, so are no 
longer entangled. The pairs that are entangled at this 
point at (3,5) and (2,6), which have become so through 
the measurement interaction. 

The exact forms of these descriptors is also interest- 
ing. We know that, because of the locality of interac- 
tions in the Deutsch-Hayden representation, the only way 
in which a descriptor can have a non-trivial dependence 
on Ha is by interacting either directly with system a, 
or with something that has interacted with it. We can 
therefore trace, not only dependencies on specific quanti- 
ties as Deutsch and Hayden did, but also trace non-trivial 
elements of a descriptor on a certain Hilbert space asso- 
ciated with another system. If we look at the descriptors 
for qubits 2 and 3 we can see that they have dependen- 
cies on many spaces, and we can trace where they have 
come from. Qubit 2 has non-trivial elements on Tii, 2,3,6- 
The dependence on Hi comes from the original Bell state 
that the pair (1,2) shared. The H3 element comes from 



'^he Bell gate interaction with qubit 3, and the He ele- 
ent from the CNOT operation with qubit 6. Similarly 
/ffe can traces the dependencies of qubit 3, with are on 
^2,3,4,5- We can also sec that qubits 1 and 4 remain 
as they were at the beginning: they have not interacted 
any further in the protocol, and this is reflected in the 
components of their descriptors. Finally we note that 
the Hilbert space dependencies of qubits 5 and 6 are the 
same as those of the qubits that they interacted with in 
the CNOT operations. 

In order to fully 'swap' the entanglement, we are going 
to need to look at the various descriptors relative to the 
four binary numbers stored in qubits 5 and 6. By con- 
struction, the pair (2,3) will be in one of the four Bell 
states, but what about the pair (1,4)? First, we need the 
four-system version of our expression (jlSp for relative de- 
scriptors. This is easily shown to be 

1l2i] = 9al2ii(l + y^(g347tm)g34nm) 



Now we use the fact that 

Wdinm) = (c^n C^m)|/3)(/3| 

where is the set of states of (3,4) relative 

to which we are finding the descriptors for (1,2). If 
mm = {|00)(00|, |01)(01|, |10)(10|, |11)(11|} then we 
have 



9l2y = 9l2-y(l ± (g3z + qiz 

i'i2tj = qi2tj{i ± (93^ - q^z 



00 

q?.zqAz) Yi 
01 

qzzqiz) Yg 



If we look at the descriptors for qubit 1 and 4, and the 
CT^-components of qubits 5 and 6, we see that the only 
nonzero components of the relative descriptors are 

qixqixq^z, qiyqiyq^izqezi qizqizq&z 

with their appropriate signs. We can therefore represent 
the relative descriptors in this case by 



qix 
qiz 



(+ H )qixqbz 

(h 1 — )qizq(iz 



which gives us the full set of descriptors relative to the 
set of measurements on (5,6) {00, 01, 10, 11}: 




qi 



If we now look at the reduced form of these descriptors 
on 7^i_4 only, we see that the simply reduced forms retain 
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all the same average values. We may therefore write the 
descriptors as 

~ay (g) CTj; q4 = (H 1 — ) -ax ^ a, 

which are the descriptors corresponding to the four 
(pure) Bell states. 

The first thing to note when analysing the flow of de- 
pendencies is that finding the relative descriptors requires 
the components q^z and qgz- Therefore, the 'swapping' 
of entanglement can only happen when these components 
are transmitted to qubits 1 and 4. That is, each qubit 
needs two bits of classical communication sent to it before 
the two qubits can become entangled. 

We therefore have the transmission of the dependencies 
of q5z,6z to qubits 1 and 4, and it is only once this has 
happened that we see the entanglement being swapped. 
Qubits 5 and 6 have dependencies on all of qubits 1 to 
4, and these are transmitted to qubits 1 and 4, possibly 
separately. The dependencies on (2,3) work out to be the 
correlation between states of (2,3) and (1,4), and those 
on (1,4) give rise to the entanglement between qubits 1 
and 4. 

Thus we see that there is nothing non-local about en- 
tanglement swapping. There is no instantaneous "trans- 
mission" of entanglement to the qubit pair (1,4) because 
of the measurement on the pair (2,3), and no superlu- 
minal signalling between qubit 1 and 4 to tell them that 



they are now entangled with each other. The entangle- 
ment dependencies (and those which give the correlation 
between the two pairs - a facet of entanglement swapping 
that is often overlooked) are transmitted to qubits 1 and 
4 separately through qubits 5 and 6. As it is only their ct^ 
components which are used, we can describe them as two 
bits of data being transmitted over a classical communi- 
cation channel, giving the result of the measurement of 
the pair (2,3) in the Bell basis. This, then, is a completely 
local description of entanglement swapping. 



IX. CONCLUSIONS 

In this paper we have developed the Deutsch-Hayden 
formalism to cover areas that it could not originally, 
specifically dealing with measurement and mixed states. 
We have also seen that, when considering the picture 
of information flow given to us by the descriptors, we 
can be confident that this is a unique picture, not de- 
pendent on the exact form chosen for those descriptors. 
Finally we have considered the entanglement swapping 
protocol, and seen that the usual description of it as cre- 
ating entanglement non-locally is incorrect, and that in 
fact the entanglement is transmitted entirely locally. We 
have therefore seen how the Deutsch-Hayden approach 
is transparent to the notions of locality and unitarity in 
quantum mechanics, in a way that the standard formal- 
ism is not. 
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